Entanglement and decoherence of N atoms and a mesoscopic field in a cavity 
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We investigate the behavior of N atoms resonantly coupled to a single electromagnetic field mode 
sustained by a high quality cavity, containing a mesoscopic coherent field. We show with a simple 
effective hamiltonian model that the strong coupling between the cavity and the atoms produces an 
atom-field entangled state, involving N + 1 nearly-coherent components slowly rotating at different 
paces in the phase plane. The periodic overlap of these components results in a complex collapse 
and revival pattern for the Rabi oscillation. We study the influence of decoherence due to the finite 
cavity quality factor. We propose a simple analytical model, based on the Monte Carlo approach 
to relaxation. We compare its predictions with exact calculations and show that these interesting 
effects could realistically be observed on a two or three atoms sample in a 15 photons field with 
circular Rydberg atoms and superconducting cavities. 

PACS numbers: 42.50.Pq, 42.50.Dv, 42.50.Ct, 03.65.Yz 



I. INTRODUCTION 

Cavity quantum electrodynamics experiments with cir- 
cular Rydberg atoms and superconducting cavities are 
well suited for the realization of tests of fundamental 
quantum processes and of simple quantum information 
processing functions 0. They make it possible, in par- 
ticular, to prepare mesoscopic quantum superpositions, 
made of coherent field components with different classi- 
cal attributes (phase and amplitude). They have opened 
the way to studies of the decoherence dynamics on these 
states, at the quantum/classical boundary Q- These 
early experiments, involving fields containing a few pho- 
tons only, were based on the dispersive atom-field interac- 
tion. The atom, off resonance with the cavity mode, be- 
haved as a state dependent transparent dielectrics mod- 
ifying transiently the cavity frequency and, hence, the 
field phase. An atom in a superposition of levels pro- 
duces then a quantum superposition of phase shifts, a 
situation reminiscent of the famous Schrodinger cat sit- 
uation. 

Much faster phase shifts can be realized through the 
resonant atom-cavity interaction. The complex Rabi os- 
cillation phenomenon in a mesoscopic field results in an 
atom-field entanglement induced by photon grainincss. 
The initially coherent cavity field is rapidly cast in a 
superposition of two components with different phases. 
This phase splitting is a mesoscopic effect that disappears 
in the classical limit of a very large field, which is then left 
unaffected by the atoms. This resonant phase splitting 
effect has been evidenced for fields containing up to a few 
tens of photons . Its coherence has been checked using 
an echo technique borrowed from NMR [5|, following a 



proposal by Morigi et al |f| . The resonant atom- field in- 
teraction thus opens the way to decoherence studies with 
large photon numbers. These experiments focused on a 
simple situation with a single atom coupled to the cavity 
mode. Recent experimental advances |g allow us to en- 
vision experiments with samples containing a well known 
number N > 1 of atoms. They would merge the concepts 
of cavity QED with the atomic ensemble manipulations 
recently put forth for quantum information processing. 
In this context, it is particularly interesting to study the 
resonant interaction of such a multi-atom sample with a 
mesoscopic field. 

In this paper, we study the resonant interaction of an 
atomic ensemble of N atoms with a cavity initially pre- 
pared in a mesoscopic coherent state. Using an appropri- 
ate mesoscopic approximation, we show that the strong 
atom/field interaction leads to an entangled atom-field 
state involving N + 1 nearly coherent field components 
with different classicalphascs, generalizing the results ob- 
tained for one atom Q- These coherent components are 
correlated with dipole atomic states, superpositions of 
the upper and lower states with equal weights. Thus, in 
the mesoscopic limit, the cavity field acts as a which-path 
detector for the atomic states interference. The periodic 
partial disentanglement of the atom-field system due to 
the transient overlap of field components is then closely 
linked to the complex pattern of quantum Rabi oscilla- 
tion collapses and revivals observed in this regime. As 
in the single atom case, early quantum revivals can be 
induced by an echo sequence, realizing a time-reversal of 
the atom- field evolution p|. The experimental observa- 
tion of these effects would shed light on the deep links 
between entanglement and complementarity. 
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This complex phase-splitting was already predicted by 
Knight and Shore 0. In the present paper, the introduc- 
tion of an effective Hamiltonian valid in the mcsoscopic 
domain enables us to capture the main results within a 
simple analytical model. This approach, originally pio- 
neered by Klimov and Chumakov 9], is also instrumental 
in the discussion of dissipation in the system. 

Dissipation in the cavity turns the entangled atoms- 
field state into a statistical mixture, destroying Rabi os- 
cillation revivals. In order to assess the experimental ac- 
cessibility of these mesoscopic quantum effects, we have 
analyzed quantitatively the influence of cavity dissipation 
on the evolution of the atom-cavity entangled state. Us- 
ing the physical insight provided by the stochastic wave 
function approach |lf| to the dissipative dynamics of the 
atoms + cavity system, an analytic formula for the de- 
coherence of the mesoscopic atoms + cavity state is de- 
rived. It generalizes to the case of TV > 1 atoms the re- 
sults previously obtained by Gea-Banacloche m the 
N = 1 case. We provide a functional expression for the 
decoherence coefficients of the entangled atoms + cavity 
state which is valid even in the presence of an echo se- 
quence used to induce an early revival of the quantum 
Rabi oscillation. The functional form of these decoher- 
ence coefficients reflects the cumulative construction of 
the imprint left by the strongly coupled atoms + cav- 
ity system in the cavity environment. These physically 
illuminating expressions can be straightforwardly gener- 
alized to compute decoherence properties during a more 
complex protocol, such as the injection of another atomic 
ensemble in the cavity shortly after the first one in order 
to probe the cavity field. 



The organization of this paper is as follows: in section 
im the model for the resonantly coupled atoms + cavity 
system is presented and its dynamics is studied in the 
absence of dissipation using a mesoscopic approximation 
in the spirit of Gea-Banacloche Q • In section IIIII dis- 
sipation of the cavity is introduced and studied analyti- 
cally using the stochastic wave function approach. Sec- 
tion IIVI presents numerical results obtained from quan- 
tum Monte-Carlo simulations. These results are used 
to discuss experimentally accessible windows for the ob- 
servation of a mcsoscopic entanglement between two or 
three atoms in a microwave high quality cavity in the near 
future. We also comment on the possibility of observ- 
ing such mesoscopic effects within the context of circuit- 
QED experiments performed with nanofabricated super- 
conducting circuits [l^. The next generation of these 
experiments will involve several qubits coupled to a cav- 
ity. Therefore, it is very natural to address the question 
of entanglement between several qubits and the resonator 
for circuit-QED devices. 



II. HAMILTONIAN EFFECTIVE DYNAMICS IN 
THE MESOSCOPIC REGIME 

A. The Tavis-Cummings model 

In this paper, the resonant interaction between N two- 
level atoms and an electromagnetic mode in a cavity is 
considered. The cavity mode is modeled by a quantum 
harmonic oscillator which, in section ITTT1 will be weakly 
coupled to an harmonic bath representing its environ- 
ment. 

Assuming that all atoms are symmetrically coupled to 
the mode, the atom-field system is conveniently described 
by the Tavis-Cummings model [l^, a spin J = N/2 
generalization of the Jaynes-Cummings model [T^. The 
interaction between the atoms and the electromagnetic 
mode is given by: 

h N 

H=f J2(Sta + SrJ) (1) 

i=l 

where Sf denote the raising and lowering operators for 
the ith atom. The energy scale associated with the in- 
teraction of one atom with the mode is hg. Because of 
the symmetric coupling, the evolution is restricted to the 
symmetric subspace, invariant under atomic permuta- 
tions, provided the initial state is also symmetric, a con- 
dition that we assume fulfilled from now on. The atomic 
degree of freedom is the spin J = N/2 representation for 
the collective su(2) generators: 

TV N 
i=l i=l 

The interaction Hamiltonian can then be rewritten in 
terms of these operators leading to the Tavis-Cummings 
(TC) model: 

H TC = ^ (j+a + J-J) . (3) 

Within this framework, the atomic ensemble behaves as 
a collective quantum object, a spin J = N/2 interact- 
ing with a quantum harmonic oscillator. A convenient 
basis in the atom + cavity Hilbert space is made up of 
tensor products of the atomic Dicke states \J,m). com- 
mon eigenstates of J 2 and J z , and the Fock states |n) for 
the harmonic oscillator. Note that the Hilbert space for 
this coupled system contains stable subspaces under time 
evolution which organize as follows: first, an infinity of 
2J + 1 dimensional subspaces H n (n > 0) generated by 
the states | J, J — I) ® |n + I) where I ranges from to 
23 . Then, a finite number of lower dimension subspaces 
indexed by — J < m < J — 1 generated by | J, m — 1} ® \l) 
where < / < J — m. 

In this paper, we focus on the mesoscopic regime in 
which the exchange of quanta between the collective state 
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of the N atoms and the cavity mode does not signifi- 
cantly alter the latter. Since the collective atomic spin 
can transfer at most N photons to the electromagnetic 
mode, this implies n 3> N, where n is the mean photon 
number in the cavity. 

B. Mesoscopic entanglement involving one atom in 
a cavity 

1. Mesoscopic approximation for the atom + cavity 
evolution 

The quantum dynamics of a single atom interacting 
with a coherent state in a cavity has been investigated 
by Gea-Banacloche Q and independently by Buzek and 
Knight ^{|. The analysis by Gea-Banacloche is based 
on the exact diagonalization of the Jaynes-Cummings 
Hamiltonian. It provides an approximate solution for the 
Schrodingcr equation with the initial condition |^> a t)<8>|o!) 
where \ipat) denotes the initial state of the two-level atom 
and | a) is a coherent state of the cavity field containing 
a mesoscopic number h = \ce\ 2 of photons (a = y/n). 

The atom + field interaction is expected to create an 
entangled state. As noticed by Knight and Shore, us- 
ing an argument based on the Schmidt theorem |8|], a 
two-level atom cannot get entangled with more than two 
orthonormal states of the field. The Gea-Banacloche 
approximate solution precisely expresses the atom+field 
state |$ (£)) at time t as a two-component entangled state. 
As discussed by Gea-Banacloche Q, this approximation 
is accurate for t <C n/g which, for n 3> 1, is large com- 
pared to the vacuum Rabi period 2n/g. It leads to: 

|tf(t)) = Ae- iatVK ' 2 \D + {t))®\i> + {t)) 

+ Be l9t ^ 2 \D.(t)) ®|iMt)> ( 4 ) 

where A and B characterize the initial atomic state (|e) in 
recent experiments 0). The atomic dipole states \D±(t)} 
are given by: 

\D ± {t)) = ^=(±e^ t /^\+) + \-)) (5) 
and the field states \ip±(t)) are: 

oo t. 

\^t))=e^e-^^m^\k). (6) 
fc=o V k - 

In the following, we will use the short hand denomination 
'Gea-Banacloche states' for these cavity states and their 
generalization to N > 1. 

2. Discussion 

Gea-Banacloche has also shown that, for times short 
compared to g^yffi. the state \ip±(t)) can be ap- 
proximated by a coherent state of parameter a±(t) = 
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FIG. 1: Schematic evolution of the entangled state for one 
atom and a mesoscopic coherent state in a cavity for real pos- 
itive a. The atomic dipole states are represented as arrows. 
The field coherent states are represented as an uncertainty 
disk at the tip of the classical amplitude. Each component 
\D±(t)) (g> \ ip±(t)) of the superposition involves an atomic po- 
larization and a field state slowly rotating in the phase plane 
at velocities ±g/4\/n~. 

e TI9 */ 4 ^a. This result is obtained by expanding y/k at 
first order in k — h around \ffi leading to: 

\ip±{t)) ~ e^ 9 *^/ 4 lae™*/ 4 ^) (7) 

in the limit t <C V^/g- Thus, the states \ip±(t)) mainly 
evolve at slow frequencies ±g/4\/fi,. We refer to JJJ) as the 
'coherent state approximation' for Gea-Banacloche states 
and, if inserted in Q , as the coherent state approxima- 
tion for the atom + cavity system. As discussed in de- 
tails in 0, this approximation breaks down for t > V^/g 
because the states \ip±(t)} undergo a slow phase spread- 
ing due to higher order terms in their expansion. They 
can no longer be considered as coherent. However, even 
if it breaks down before the mesoscopic approximation, 
the coherent state approximation provides a nice intu- 
itive and pictorial support for visualizing the system's 
evolution. 

With this image in mind, it is useful to draw on the 
same diagram the motion of the average atomic polariza- 
tion d±(t) = (D±(t)\ a \D±{t)) in the equatorial plane 
of the Bloch sphere and the motion of a*±{t) in the 
Fresnel plane. The result is depicted on figure ^ for 
ip = 0: the corresponding vectors rotate at angular ve- 
locity ±g /4a/^, small compared to the classical Rabi fre- 
quency g\ffi. 

In the limit n — > oo and gt -C 1, both states \ip± (t)} are 
close to | a) for t <C g~ x , meaning that the cavity mode 
is barely affected by the atoms. In this regime, the cav- 
ity state factors out and the atomic polarizations \D±(t)) 
coincide with the atomic spin-1/2 eigenstates along the 
x direction. These atomic states interfere resulting in 



the classical Rabi oscillation phenomenon. Rcmembt 
that the period of classical Rabi oscillations is of the o: 
der g~ x j\ffi. Therefore, in the classical limit n — ► +cx 
more and more oscillations take place before the motic 
of Gea-Banacloche states in the phase plane has any me; 
surable consequence. 

In the mesoscopic regime (fixed n ^> 1), the stal 
of the cavity is altered by the atom. Interferences bi 
tween atomic polarizations \D±(t)) can only be observe 
when \ijj + (t)) and \ip-(t)) overlap. As explained abov 
at very short times, these states are still close to tl 
initial coherent state \a). When the phase separatic 
between \tp + (t)} and \ip-{t)) due to their slow rotatic 
in phase space is larger than the quantum phase flucti 
ations in these coherent components (gt/2y/E ~ 
the cavity field behaves as a bona fide "path detectoi 
for the atomic polarizations and Rabi oscillations disaj 
pear. The Rabi oscillation collapses after a time of the 
order of the vacuum Rabi oscillation, after \ffi classical 
oscillations. 

The Rabi oscillation signal reappears when \tfj + (t)) and 
\ip-(t)) overlap again. This happens for gt/2y/E ~ 2ir. 
During this overlap, the disentanglement of the atom + 
cavity state erases the information stored in the cavity 
about the path followed by the atomic degrees of freedom. 
This 'quantum eraser situation' leads to a revival of Rabi 
oscillations. Rabi oscillation revivals in the mesoscopic 
regime are thus a direct application of the complemen- 
tarity concept ^(|. Figure |3 shows, as a function of the 
dimensionless time = <)tf/2\/n, the first spontaneous re- 
vival of the Rabi oscillation signal obtained by numerical 
integration of the Schrodingcr equation for one atom ini- 
tially in the excited state and coherent states of 14 and 
40 photons in average. 



3. The echo protocol 

The echo protocol proposed by Morigi [l| aims at test- 
ing the coherence of the atom + cavity state by a time 
reversal operation. A percussional echo pulse is applied 
to the atom at time t^. It corresponds to the unitary 
operator U n = ia z . The evolution for the atom + cavity 
system up to time t > is then given by 



U(t) 



-j(t- 



where H is the Jaynes-Cummings Hamiltonian. 
Ul = 1 and U„HU n = -H, we get: 



U(t) = U*. 



-i(t„-t)H/h 



rH/fr 



(8) 
Using 



(9) 



Therefore, right after the echo pulse, the Gea-Banacloche 
states reverse their evolution and recombine at 2t T lead- 
ing to an induced Rabi oscillation revival. This induced 
revival may occur at much shorter time than the 'sponta- 
neous' revival, making its experimental observation much 
easier, as shown recently £|. Moreover, in the absence of 



(a) 




(b) 




FIG. 2: Spontaneous revivals of the Rabi oscillation signal 
for one atom with initial condition (|e) = |m = 1/2)) ® | a) 
and a coherent state with mean photon number (a) n = 15 
and (b) n = 40. The solid line shows the probability P{(j>) 
for finding the atom in |e), as a function of the dimensionless 
time 4> = gt/2vft, computed using numerical integration. The 
dashed curves show the upper and lower envelopes predicted 
by our mesoscopic approximation. 



decoherence, the induced revival should occur with unit 
contrast. The influence of decoherence could thus be, in 
principle, directly assessed from the measurement of the 
induced revival contrast. 



4- Towards atomic ensembles 

In this paper, we are interested in studying the reso- 
nant interaction of an atomic ensemble containing N > 1 
atoms with a mesoscopic field in a cavity. Invoking again 
Schmidt theorem [j| , we expect this resonant interaction 
to create an entangled state with 2 J + 1 = N + 1 or- 
thonormal components. As in the one-atom case, partial 
disentanglement of this state will lead to spontaneous re- 
vivals of Rabi oscillations. 

An analytic diagonalization of the Tavis-Cummings 
Hamiltonian can be obtained for N = 1,2,3 but not for 
greater values of N. Moreover, as will be clear from forth- 
coming sections, the analytical diagonalization for these 
values of N does not enlighten the dynamics of the sys- 
tem. In particular, for N > 1, the explicit expressions of 
exact eigenstates in the resonant Tavis-Cummings Hamil- 
tonian depend on N. This direct approach thus cannot 
be used as a convenient starting point for an approximate 
solution of the Schrodingcr equation for N > 1. 

Our approach, dcvcloppcd in the next section, relies 
on an effective Hamiltonian which, in the mesoscopic do- 
main, provides an excellent approximation to the Tavis- 
Cummings Hamiltonian. It provides a unified vision of 
the dynamics for all values of N and as such, it is a 
good starting point for analyzing the dynamics in the 
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mesoscopic domain. As we shall see, in this framework, 
the dynamics of the resonant Tavis-Cummings model can 
then be described in the spirit of the Gea-Banacloche ap- 
proach. 



C. Effective dynamics in the mesoscopic domain 

1. Effective Hamiltonian 

First of all, let us remark that any initial state of 
the form \J,mo) ® \a) in the mesoscopic domain mainly 
spreads over (2J + l)-dimensional stable subspaces Ti n 
for values of n around n. The core of our approach is to 
replace the Tavis-Cummings Hamiltonian © acting on 
subspaces Tl n by an effective Hamiltonian in which the 
n dependence factors out. It appears that the main n 
dependence of J3J scales as y/ri for large values of n. As 
in the N = 1 case, this non-linearity leads to the collapse 
of the Rabi oscillations and the discrete character of the 
spectrum leads to spontaneous revivals. 

In order to describe our ansatz for the effective Hamil- 
tonian, it is convenient to remark that each subspace 
Tin can be turned into a spin- J representation of su(2). 
Let us introduce new operators J^ 1 and J z . With 
the notation \Zm ) = \J,m) ® \n + J — m), these new 
operators simply act on these states in the same way 
as standard su(2) generators act on the \J, m) states: 
J^Z^) = y /J{J+\)-m{m±l)\Z%l 1 ) (see fig. 0. 
The operators a J + and J~ then act on the states 
\zt ] ) as: 

a.J+\Z^) = y/n+J-mJ+\Z$) (10) 
<JJ-\Z^) = y/ n + J-m + lJ-\zM). (11) 

We then note that computing the evolution of a state 
| J, m) ®\a) in the mesoscopic regime requires considering 
values of n close to n > iV. The variation of \Jn + k 
for < k < N + 1 is small for n > N (of t he order of 
N/y/n). We thus drop the m-dependence of \Jn + J — m 
and vn + J — m + 1 by replacing them by yjn + c, where 
0<c<A^+lisa constant to be discussed in the next 
paragraph. This leads to an effective Hamiltonian of the 
form: 



#eff - y ( J + +J~) = hgV^+~c J x . (12) 

This Hamiltonian, already derived by Klimov and Chu- 
makov and used to study the squeezing of light by 
an atomic ensemble jl7j . shares some features with the 
expected classical dynamics, driven by an effective field 
along x. Here however, photon emission and absorption 
are taken into account through the fact that J ± changes 
the photon number (see fig. 0J. At fixed n, the eigenval- 
ues of this effective Hamiltonian are equally spaced, as 
predicted by [T^ | for the Tavis-Cummings Hamiltonian in 




-3/2 -1/2 1/2 3/2 

FIG. 3: Schematic view of the Hilbert space for N = 3 
atoms. The oblique grey lines represent the stable subspaces. 
The action of J± on Hs is depicted. 



the large n limit. All the n dependence of this effective 
Hamiltonian is contained in the \Jn + c factor. 

Of course, there is an ambiguity in the choice of < 
c < N + 1 but we shall see that (i) for N — 1 choosing 
c = 1 reproduces the results of section III Bl and (ii) for 
N > 1, changing c only affects the rapidly oscillating 
part of the Rabi oscillation signal. It does not change its 
envelope which is precisely the information we hope to 
extract from the effective Hamiltonian. 



2. The mesoscopic approximation 

Using this effective Hamiltonian, it is possible to study 
the evolution of a state = \J, m )x ® \oi) where 

J x \J,m) x = m \J,m) x . An approximate solution for the 
Schrodinger equation shows that this state remains fac- 
torized (see appendix lAl for details): 

!*£(*)> = e~ im9tH \D m (t)) ® \i>m{t)) (13) 

where the state of the electromagnetic mode is of the 
form: 

|V> m (t)) = e im3tVW e-"/ 2 V ^= e-^^k) , (14) 



which we call, as above, a Gea-Banacloche state |7j • The 
atomic polarizations generalize the ones found by Gea- 
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Banacloche in the spin 1 /2 case: 
,/ 

\D m (t))= J2 e-'o^-^'y^iR-^m^l^m'). 

m / — — J 

(15) 

where R denotes the rotation matrix R m ,m' = 
(J,m'\e l7rJV/2 \J,m). Note the presence of the classical 
Rabi frequencies ?n<?|a| corresponding to the quantum 
beat between spin eigenstates along the x direction of 
the effective classical field. The average angular mo- 
mentum d m (t) = (D m (t)\ J \D m (t)} slowly rotates in the 
equatorial plane of the Bloch sphere at angular veloc- 
ity gm/2\/E. The parameter c appears in these atomic 
polarizations only and, for N = 1, the Gea-Banacloche 
results are exactly recovered for c = 1. 

Starting from state |vP(0)) = \J, m) g) \a), an entangled 
state with N + 1 components is obtained: 

J 

|*(t))= Rmo,rne- imgtlal \D m (t))®\lP m (t)}. (16) 

m— — J 

As in the one-atom case, the entangled state <|16ll can be 
viewed as the result of the ideal measurement of the spin 
by the mesoscopic field in the cavity. 

3. The coherent state approximation 

As in the single-atom case, the state |^ m (f)) can be 
approximated by a coherent state of complex amplitude 
ctmit) = e _lm9 */ 2v ^a. This approximation holds in 
the limit t <C (</|m|) _1 -\/^- At longer times, typically 
y/E/g\m\, the field state gets deformed as the \ip±) states 
in the N = 1 case. 

With this image in mind, it is useful to draw on the 
same diagram the motion of the average atomic polar- 
ization d m (t) in the equatorial plane of the Bloch sphere 
and the motion of a m (t) in the Fresnel plane, generaliz- 
ing the phase space representation used above. The main 
difference with the N = 1 case is the appearance of N+ 1 
frequencies and field states instead of two. Here also, the 
phase of the coherent state plays the role of a pointer 
measuring the angular momentum of the collective spin 
along the x direction in the n 3> 1 limit. Larger angular 
momenta lead to larger angular velocities. As we shall 
see now, this complex atoms + cavity entangled state 
leads to a rich pattern of spontaneous revivals of Rabi 
oscillations. 



D. Partial revivals of Rabi oscillations 

1. General picture 

Rabi oscillations of the atomic populations provide a 
nice way to probe the degree of entanglement of the atom 
+ cavity state. 



In the classical limit (n — > oo), the electromagnetic 
field state factors out and quantum interferences between 
the various atomic polarizations \D m {t)) can be observed. 
They are the Rabi oscillations for the quantum spin J in 
a transverse classical field. 

In the mesoscopic limit, the electromagnetic mode is 
altered by the atom. Interferences between atomic po- 
larizations \D m (t)) can only be observed when the cor- 
responding field states \ip m (t)) overlap. At very short 
times, the various components \tp m (t)) are still close to 
the initial coherent state \a) and Rabi oscillations show 
up. When the various Gea-Banacloche states split apart, 
the electromagnetic field becomes a good "path detec- 
tor" for the atomic polarizations and the Rabi oscilla- 
tions collapse. This is again a complementarity effect, 
the field storing a which-path information about the in- 
terfering atomic states. Rabi oscillations reappear when 
this which-path information is, at least partially, erased, 
i.e. when some of the Gea-Banacloche states overlap 
again. In the N > 1 case, the atoms + cavity state is a 
superposition of N + 1 factorized components rotating at 
different velocities and a rich spontaneous revival pattern 
is expected. 

The angular velocity of the Gea-Banacloche states 
\ipm(t)) in the Fresnel plane is here —gmj2\ffi 1 suggest- 
ing to associate with each Rabi oscillation revival a non- 
empty subset £ of {1, . . . , 2 J} which is the list of abso- 
lute values of differences of the indices m of those Gea- 
Banacloche states that overlap during the revival un- 
der consideration. Given such a subset £ , the revival 
is built from contributions of pairs of states \tp m+ (t)) 
and \tfj m _(t)) such that |m+ — m_| G £. Note that, 
when a pair \ip m+ (t)) and \ipm-(t)} overlaps, all pairs 
that have a positive or negative integer multiple value of 
m_ — m + also overlap. Therefore if q belongs to such 
a subset, its multiples also do. This is the only con- 
straint on the subsets £. The Rabi oscillation revivals 
are therefore classified by the greatest common divisor 
gcd(£) of the elements of £. The first time of occurrence 
of the spontaneous revival associated with £ is t such that 
gt/2y/E = 2-7r/gcd(£). Replicas of this revival will occur 
at integer multiples of this fundamental time. Note that 
there are N + 1 — gcd(£ ) pairs of Gea-Banacloche states 
that verify \m + — m_| = q. Remember that <fi = gtjl^ffi 
is the dimensionless time which characterizes the slow 
motion of Gea-Banacloche states. In general, the con- 
trast of replicas will be reduced because of the spreading 
of the Gea-Banacloche state (especially if they occur af- 
ter gt/2y/n > 2tt). For all values of N, the set {1, ... , N} 
corresponds to a complete revival involving the recombi- 
nation of all Gea-Banacloche states at <fi = 2ir. 

As an example, let us consider the case of three atoms. 
The corresponding Rabi revival patterns are depicted on 
Figure ^] The first revival is obtained when m = ±3/2 
states overlap for <fi = 2tt/3 (associated subset £ = {3}). 
It is partial (contrast is lower than one) since only two 
atomic polarizations take part in it. The next revival 
appears for (f> — n when m = 3/2 & m = —1/2 and m = 
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Im(a) 




Im(ct) 



Im(a) 



Re(a) 



Re(a) 



Im(a) 



Re (a) 



Re(a) 



Within our effective Hamiltonian approximation, the 
time dependence of the atomic polarization factor 
VnT-Jml is a phase. Modulation factors for the revivals 
come from the overlaps of cavity mode states: 

Tl m+ , m _ (t) = (lp m - {t)\1>m+ (*)) . (20) 

This expression only depends on q = m+ — TO- . Finally, 
the Rabi oscillation signal is 

2J 

P m {t)= K g (t)4'"°' m) e" l ^e^ Mtv/ * (21) 

q=~2J 

where 



m_ — m+ — q 



FIG. 4: Position of the Gea-Banacloche states at the times 
of first occurrence of spontaneous Rabi oscillation revivals for 
N — 3 (J = 3/2): (a) Collapse of the Rabi oscillations when 
the various components split apart (b) First spontaneous re- 
vival for <f> = 27r/3 (c) Second spontaneous revival for <f> — tt 
(d) Complete spontaneous revival involving all atomic polar- 
izations for cj> — 2tt. 

—3/2 & m = 1/2 overlap separately (£ = {2}). At = 
47r/3, m = 3/2 and m = —3/2 recombine again leading 
to a partial revival which is a replica of the first one {£ = 
{3}). Finally the complete revival involving quantum 
interferences between all four atomic polarizations takes 
place at <fr = 2ir {£ = {1,2,3}). Note that the first 
example of a non-trivial revival involving several slow 
frequencies before the complete revival occurs for N = 4 
at <jy = tt (S = {2,4}). 

2. Rabi oscillation envelopes 

Preparing J + too atoms in the excited state, the prob- 
ability of detecting J + m ones in the excited state and 
J — to in the ground state is given by (q = m + — m_ ) : 

Pm(t)= £ e-^^V^l{t)n m+ . m ^ (17) 

where m± run from — J to J and 

= Rmo,m + R* mo , m _ (J, m \D m+ (t))(D m _(t)\J im ) 

(18) 

contains the matrix elements of atomic polarizations. 
The scalar products (J, m\D m± (t)) can also be expressed 
in terms of the rotation matrices (see appendix f° r 
explicit expressions): 

(J,m\D m± (t)) = (ir% ± , m e-^ c - J+m ) m ± t / 2 ^. 

(19) 



This expression separates the rapid frequencies gq\fn 
from the mesoscopic slow frequencies gq/2VE. The Rabi 
signal then consists into a rapidly oscillating signal slowly 
modulated in amplitude and phase. Expression l|21() can 
be used to find an approximate analytic expression for 
the upper and lower envelopes of the signal. 

Let us illustrate this point on the signal obtained for 
mo = m = J (denoting A q J,J ^ — A q ) which is plotted in 
the forthcoming figures. 

Outside the spontaneous revivals, the contribution of 
the q ^ terms in l|21(l vanishes. The base line of the 
Rabi oscillation signal is thus Aq. Obtaining the en- 
velopes is trivial for N = 1 since there is exactly one 
value of q involved: q r = 1. The slowly varying phase 
e -i9<?tc/2\/fi s i m ply shifts the rapid oscillation without 
changing its envelope. This analysis is correct also for 
N > 1 in the case of a revival involving exactly one 
frequency corresponding to q r £ {1,...N} (q r = 1 for 
N = 1). For these revivals only, the upper and lower 
envelopes P+ and P_ take the form: 

P ± (t)=A ±\n qr (t)A qr \. (22) 

Thus, simple revivals are symmetric with respect to 
the flat signal Aq and involve only one rapid frequency 
gq r V^- 

The analysis turns out to be more involved when sev- 
eral frequencies are involved. A first example of this situ- 
ation is the initial collapse of Rabi oscillation (i < 2jr/g) 
for N > 1. Nevertheless, the envelopes can be obtained 
using the exact expression for classical Rabi oscillation 
in a field of amplitude \/Tl: in this limit, the probabil- 
ity for detecting all atoms in the excited state is given 
by P c (t) = cos 2N {gtVn/2). Its maxima occur at times 
2i\kj g\ffi for integer values of k and its minima occur for 
half integer values of k. Substituting these values in l|21(l 
for large values of fi provides the values of the rapidly 
oscillating term to be used to fit the maxima (upper en- 



8 



vclopc) and the minima (lower envelope). This leads to: 
P+(t) = A + J2\K q (t)A q \ (23) 

q^O 

P_(t) = A + J2(-V q \K q (t)A q \. (24) 

Note that for N > 1, the envelope is not symmetric with 
respect to the flat signal Aq. 

Let us now turn to the complete revival which takes 
places around tn = Ait^/fi/ g. Near this revival, the Rabi 
oscillation signal takes the form (t = tn + t) : 

P(t) = A q K q (t R +T) e~ 1 ^ e -im{rVfi+2<c+2n)) 

q^Q 

(25) 

We first note that the rapidly oscillating phases are 
shifted in time by 2ir(c + 2n)/\ffi. This time shift does 
not affect the low frequency modulating terms A q lZ q {t). 
Within the coherent state approximation, the overlap fac- 
tors Hq(tu+T) can be approximated by 1 for |t| < 2w/g. 
This means that, in the classical limit and within the co- 
herent state approximation, close to the complete revival, 
the Rabi oscillation signal has the same fast oscillations 
than near t = 0. This suggests to use (|23|) and (|24ll as 
upper and lower envelopes. Because the Gea-Banacloche 
states are getting deformed over a time scale ta, these 
expressions only provide an approximation to the real 
envelopes of the theoretical signal (|25|) . This approxi- 
mation assumes that the overlap factors lZ q (tji) do not 
depend on q. Because \TZ q (tn) — 1| goes to zero as 1/n 
in the large n limit, the accuracy of l|23|) and l|24|l as ap- 
proximate upper and lower envelopes for the main spon- 
taneous resurgence increases with increasing n. 

To summarize, <|23|) describes the upper envelope of all 
revivals. The lower envelope is described by P—(t) in eq. 
(|22|l for revivals involving only one frequency such as the 
ones occuring at gt/2\ffi = 2ir/q where [N/2] < q < N 
and the lower envelope of the main revival which occurs 
at gt/2\ffi = 2tt is described by (|24|l . In practice, only 
the lowest values of N (N < 3) may be easily reachable in 
the Rydberg atoms experiments and therefore, the only 
revival involving more than one value of |m + — m_| is 
the complete one. 



3. Numerical results 

All numerical results in this paper are presented in 
terms of the dimensionless time <fi = gt/2\ffi associated 
with the slow evolution induced by the atom + field in- 
teraction. 

Let start by considering the case of N = 1 atom. 
Figure shows the comparison between the analytic 
envelopes (122ft and an exact numerical solution of the 
Schrodingcr equation for different values of h. As ex- 
pected, the mesoscopic approximation becomes better 



and better as n increases. In this case, the upper and 
lower envelopes of the mesoscopic approximation signal 
are obtained by setting q r = 1 in (|2*2)l . 

Let us now consider the case of three atoms. Figure 
presents a comparison between the results of an analytic 
exact diagonalization of the Tavis-Cummings Hamilto- 
nian and the mesoscopic approximation. Correspondence 
with the revivals described in fig. 0]is indicated. The up- 
per and lower envelopes l|23[) and (|24() are depicted. Fig. 
[SJshows that although the effective Hamiltonian does not 
fully reproduce the exact signal, it does reproduce the 
amplitude and the positions of the revivals in a satisfac- 
tory way. 

As expected, eqs. <{2*3")l and (|2*4*|) effectively describe the 
upper and lower envelopes of the signal during the early 
collapse of Rabi oscillations. They fit also rather well 
with the first complete revival (cf> ~ 2ir). But they fail 
for the revival at 4> ~ (|23|) corresponds to the upper 
envelope but (|24|) does not. This is not surprising since 
this revival is due to q = 2. In this case, formulas 12- il) 
and (|24|l weigh the q = 2 contribution with the same 
sign. Eq. 12 211 with q r = 2 would be more appropriate to 
describe the envelopes near this revival. 

Finally, the value Aq of the probability between re- 
vivals obtained from the effective Hamiltonian differs 
from the one obtained from the numerical solution of the 
Schrodingcr equation in the case N = 3. This is a finite n 
effect arising from the choice of the effective Hamiltonian 
(|12H . It can be checked that this difference vanishes as 
?i -1 / 2 in the classical limit n — * +oo. 

Comparing Rabi oscillation revival patterns at fixed n 
for various N shows that using two or three atoms instead 
of one induces an earlier spontaneous revival because the 
extreme Gea-Banacloche states (m — ±J) move faster 
than the ones associated with m = ±1/2. But the weight 
of Rabi oscillations generated by high \m\ polarizations 
quickly decreases with N. Fig. suggest that the first 
spontaneous Rabi revival for N = 2 and the second one 
for N = 3 could be good candidates for the observation 
of spontaneous Rabi revivals. Of course, dissipation in 
the cavity leads to smaller Rabi oscillations as we shall 
see in the next section. 



III. DISSIPATIVE DYNAMICS 

A. Stochastic wave function approach to quantum 
dynamics 

1. General principle 

Within the context of cQED experiments performed 
with Rydberg atoms in microwave cavities, dissipation 
almost exclusively originates in cavity losses. They are 
extremely low since the quality factor Q of the cavity is of 
the order or higher than 10 8 . Dissipation can be modeled 
through the coupling of the cavity mode to an harmonic 



9 




valid down to T 
the form: 



K and, in the present case, takes 



FIG. 5: Spontaneous revivals of the Rabi oscillation signal 
P = Pm-.j for TV = 3 atoms (J = 3/2), as a function of 
the dimensionless time <j> = gtjlyffi. The initial state |m = 
3/2} (g> \a) with an average photon number n = \a\ 2 = 15. 
The signal is computed using the mesoscopic approximation 
(grey line), its upper and lower envelopes (plain lines) and 
the exact solution (plain line with rapid oscillations) in the 
dissipationless case. Letters (a), (b), (c) and (d) refer to the 
overlaps of cavity field states for each revival depicted in fig. 
0| As expected, spontaneous revivals (b) and (c) are symmet- 




ry 3 i itt i 
— = -- [H,p\+^a.p.a 



t_Z 



— {a) a. p + p. a) a) (26) 



where H denotes the Tavis-Cummings Hamiltonian 
Note that switching to an interaction representation for 
the atoms and the cavity mode does not modify the form 
of the dissipative terms and simplifies the Hamiltonian 
part. In principle, cq. (|26|l can be solved numerically 
in order to obtain the quantum dynamics. However, an 
analytical ansatz for the reduced density matrix can be 
found within the mesoscopic approximation. As we shall 
see in the next section, this ansatz is conveniently de- 
rived using an alternative but equivalent approach to the 
dissipative dynamics of the atoms + cavity system: the 
quantum jump approach [Toj . 

The basic idea underlying this approach is to consider 
that the environment of the system is continuously mon- 
itored so that any emission or absorption of quanta by 
the system can be assigned a precise date. Each time 
such an event occurs, the system undergoes a quantum 
jump. Between these jumps, its evolution is described by 
an effective Hamiltonian that describes both its intrinsic 
dynamics and the acquisition of information arising from 
the fact that no quanta has been detected. The proba- 
bility rates for the various quantum jumps are directly 
obtained as averages of in the state under consid- 
eration where the Li denote the quantum jump operator 
(here, only L = y^a is present at zero temperature). 
The reduced density matrix is then recovered by averag- 
ing over the set of stochastic trajectories associated with 
a large set of quantum jumps sequences. The weight of a 
given trajectory can be directly related to the dates and 
types of the various quantum jumps. 

This method proves to be very convenient numerically 
since the number of variables involved is of the order of 
the dimension d of the system's Hilbert space whereas 
it scales as d 2 in the master equation approach. Note 
that the quantum jump approach is the only adequate 
formalism for studying the behavior of a single realization 
of the quantum system. 



FIG. 6: Spontaneous revivals of the Rabi oscillation signal 
for n = 15, without dissipation, computed from a numerical 
integration of the Schrodinger equation for (a) TV = 1 atom, 
(b) N — 2 atoms, (c) N = 3 atoms as a function of the 
dimensionless time 4> = gt/2\/n~- The initial condition is m = 
N/2 (all atoms excited). 



bath with very short memory. Within this framework, 
the dynamics of the coupled atom + cavity system can 
be described by a master equation for its reduced density 
matrix. The master equation is valid over time scales 
much larger than the memory time r c of the bath. In the 
weak dissipation limit, which is realized here, it is still 



2. Decoherence of coherent state Schrodinger cats 

Before applying this method to our problem, it is in- 
structive to recall how the dissipative dynamics of an 
harmonic oscillator can be described in this way. In par- 
ticular, we shall review how the decoherence scenario for 
a superposition of two coherent states can be recovered 
within this framework since it will prove to be useful in 
our cQED problem. At T = K, the system can only 
emit quanta. The stochastic dynamics of the quantum 
state is then described as follows: 

• During a small time interval r 3> r c , the probabil- 
ity for a quantum jump is p T = 'jT(ip(t)\<ra\ip(t)) 
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and the state after such a jump is \ip(t + )) = 
y/rra\tp(t))/^/p^. 

• Between jumps, the effective non-hermitian Hamil- 
tonian is given by 



ft, H g = Ll/qq T a — i—a T a 



(27) 



and between t and t+dt the state evolves according 
to 



\ip{t + dt)) 



(l - ih^H^dt) \>>p(t)) 

y/Tl (t,dt) 



(28) 



where the probability Tlo(t,dt) that no quantum 
jump occurs between t and t + dt is given, in the 
present situation, by Ho(t, dt) = 1 — pat- 

The evolution of a single coherent state \a) is well known: 
apart from a well defined global phase factor, a quantum 
jump docs not affect it. A coherent state remains coher- 
ent during the evolution between quantum jumps but its 
amplitude decreases exponentially due to the acquisition 
of information by negative measurements (no photon es- 
capes) between the jumps: a(t) = ae M °*~ 7 */ 2 . 

Let us now consider a superposition of two distinct 
coherent states with the same average number of quanta 
n but with a relative phase 9 in phase space: 



Thus, the decoherence coefficient is the characteristic 
function for the probability distribution of the number 
of quantum jumps N[Q,t] between and t. The oscilla- 
tor reduced density matrix at time t is given by: 

p(t) = \a(t))(a(t)\ + \a(t)e ie )(a(t)e ie \ (32) 
+ V{t) \a{t)e 10 ) (a{t)\ + V(t)* \a(t)) (a{t)e l9 \ 

where T>(t) denotes the average of e^ ''! 6 * over all 
stochastic trajectories. When yt > 1, the reduction of 
the components amplitude has to be taken into account 
in order to get the probability distribution for a sequence 
of quantum jumps dates < t± < . . . t p < t. This step is 
necessary to recover the full solution of the master equa- 
tion (|26() . These computations are recalled in appendix 
IO But for t <C 7 _1 , relaxation of energy has almost not 
occured and we can assume that the average number of 
quanta in the two coherent components of the superpo- 
sition is still equal to n. Therefore, emission of quanta 
is a renewal process with a distribution of waiting times 
given by ip(r) = fvy e _ ™ 7T . Decoherence by a sequence of 
quantum jumps obeying a renewal process has been re- 
cently studied in full generality |l9j . In the present case, 
since -/V[0, t] is distributed according to Poisson law with 
mean value fi^/t, the result is given by: 



■ D ( t ) = (e iN ^ e )=exp(n 7 t(e i9 - 1)) 



(33) 



Wo-)) 



i 



V2 



(\a) + \e w a 



(29) 



The jumps probability during time r is then given by 
p T = n^T (l + 2$l((a\ae ie ) e te )) which simplifies to n-/T 
as soon as \a) and \a e l6 ) are well separated so that their 
overlap can be safely neglected. Under this assumption, 
the state after a jump is given by: 



and leads to the same predictions as the direct solution 
of the master equation. 



B. Decoherence in the mesoscopic approximation 



|V(0+)) = -^(e^)|a) + e 



In the present unfolding of the master equation, each 
quantum jump introduces a phase factor e in the quan- 
tum superposition whereas each component remains a 
coherent state with the same parameter. Decoherence 
occurs because the number of jumps in a given time in- 
terval varies from one stochastic trajectory to the other. 
This has already been noticed in [l£| using an unfold- 
ing of the master equation based on a continuous time 
measurement through an homodyne detection of the field 
leaving the cavity. The present scheme leads to the same 
final results but is more suited for our study of dissipation 
on the atoms + cavity dynamics. 

Denoting by {t\, . . . , t p } the successive dates of quan- 
tum jumps (0 < t% < ... < t p < t), the final state 
associated with this sequence of jumps is given by: 



Let us now turn to the dissipativc dynamics in the 
atoms + cavity problem. Because Rydberg atom exper- 
iments are performed over time scales rather short com- 
;(Arg(a)+e) i ae i9\ \ f%Q\ pared to the dissipation time (jt < 0.1), we shall look 
/ at the dissipative dynamics at short times when energy 

dissipation can be neglected (-ft <C 1). In order to make 
an explicit connection with the work of Gea-Banaclochc 
who has studied the effect of dissipation for one atom at 
arbitrary times the case of longer times (7^ > 1) 
is discussed in appendix [Dj Let us finally mention that 
the case of an atomic ensemble has also been considered 
within the framework of master equation |20T | which, in 
our opinion, does not clarify the dissipative quantum dy- 
namics of the cavity + atoms system as much as the 
quantum stochastic trajectories method discussed below. 



1. Evolution along a single stochastic trajectory 

Inspired by the dissipationless case, we will focus on 
the evolution of factorized states of the form 



\i>{t lt ...,t p }(t)) 



1 

71 



ae 



-7*/2\ 



'>■) 



(31) 



!*£)= E KkJJ^") 



(34) 



m" — — J 
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which, within the mcsoscopic approximation, remain fac- 
torized in the absence of dissipation. Strictly speaking, 
as noticed by Gea-Banacloche for TV = 1 and as proved in 
appendixiDl this is not true in the presence of dissipation. 
However, in the limit -C 1, the dissipative dynamics 
can still be formulated in terms of factorized states. Let 
us sketch the argument that justifies this assertion. We 
refer the reader to appendix ID1 for details. 

Between quantum jumps, the state of the atoms + 
cavity system evolves according to (|28l) using the non- 
hcrmitian Hamiltonian l|27(l which takes into account the 
information gained by observing that no quanta is emit- 
ted between two jumps. Because the atoms and the cav- 
ity mode are coupled, the whole atoms + cavity state 
should be affected by this information gain. But, in the 
present strong coupling situation (g 3> 7), we expect 
the atoms to be mainly driven by the cavity and not 
by this indirect information gain. Next, dissipation in- 
duces an exponential decay of the average photon number 
n{t) = n e~ 7 * while keeping the photon number distribu- 
tion Poissonian. We shall thus neglect the decay of n{t) 
for jt <C 1. The resulting evolution for the atoms + 
cavity case is the same as in the dissipationlcss case: 

|*£(t)) ~ e -teWS \D m (t)) ® |iMt)> . (35) 

Let us now discuss the effect of a quantum jump on 
this state. Contrarily to coherent states, each Gea- 
Banacloche state \ip m (t)} does not remain invariant under 
the action of a quantum jump operator since: 

00 k 

a\i> m (t)) = e - fi /2^ a |_ e -^v^+T| fc) (36) 

fc=0 

00 fc 
= e~"/ 2 a e « Tt ff*(v / fc-Vfc+l) ^_ e -igmtVk n\ 

^ fc! 

k=0 

The phase factor cxp [imgt(s/k + 1 — Vk)] a priori de- 
pends upon k. But expanding y/k + 1 — vfc in 
powers of (fc — n)j\ffi shows that, at first order, 
exp (imgt(y/k + 1 — Vk)) is indeed independent of fc. Us- 
ing this approximation, the action of the annihilation op- 
erator reduces to the multiplication by a phase: 

a |VU*)) - a e-"^/2v^ W) . (37 ) 

It can be shown that an expansion to the next order dif- 
fers from this expression by 0((mgt/fi) 2 ). Thus, l|37|) 
can be considered valid approximation in the do- 
main t <C V"/<7 which is in the domain of validity of the 
mcsoscopic approximation (f <C n/g). 

Let us now consider the evolution of state (|34(l along a 
single stochastic trajectory. Just before its first quan- 
tum jump, provided it happens at time t\ such that 
n(ti) ^ N, we still have a factorized state of the form 
(|35J) . The effect of a quantum jump occuring a time t\ is 
to extract a phase e %6m<ytl ^ where O m (ti) is the argument 
of the amplitude of the quasi-coherent state IVVn (All- 



iterating this argument shows that, in a stochastic tra- 
jectory with quantum jumps occuring at times < t\ < 
...< tp <t, remains factorized but gets an extra 

phase 6 m (t\, . . . , t p ) associated with the quantum jumps: 

t p }; t) = e ie ^ *»> e-«^|I? m (i)>®|V> ro (*)) 

(38) 

and 6 m (ti, . . . , t p ) = J2j Qm(tj)- Exactly as for the case 
of a superposition of coherent states of an harmonic os- 
cillator analyzed above, the accumulation of random rel- 
ative phases in front of the Gea-Banacloche states leads 
to the decoherence of the Schrodinger cat state created 
by the strong atoms + cavity coupling and to the disap- 
pearance of Rabi oscillations. 

Before discussing the average over all stochastic tra- 
jectories, it is worth mentioning that the present discus- 
sion remains valid even in the presence of the echo pulses 
introduced in section iBBl A 7r-pulse instantaneously re- 
verses the dynamics of the atoms + cavity system. After 
a single 7r-pulse at time t w , the Gea-Banacloche states in- 
vert their motion and start refocusing. The deterministic 
evolution of the atoms + cavity system is then described 
by a time-reversed evolution of the dissipationless mo- 
tion. Therefore, the effect of any subsequent quantum 
jump is still to extract a phase corresponding to the po- 
sition of the quasi-coherent Gea-Banacloche state in the 
Fresnel plane at the jump time. 



2. Average over stochastic trajectories 

To deal with all these situations at once, let us denote 
by m (t) the time-dependent phase of \tjjm(t)}, not as- 
suming any particular form. The decoherence coefficient 
for the two states (t)) considered here is thus given 
by the average over sequences of quantum jumps: 

^ + ,»™j4 ,z,(A9)(t,) ) ( 3Q ) 

where (A8) = (6 m+ — m _)(t) and the t\ are the dates of 
the successive quantum jumps occuring between and t. 
This coefficient now depends in a functional way on the 
two trajectories 1 1— > 9 m± (t) and should be called a deco- 
herence functional in reference to the work of Feynman 
and Vernon [2l]]. Its definition (|39|l generalizes (|33fl to 
the case of a time dependent Ad. Since we assumed that 
jt <C 1, the statistics of waiting times between quan- 
tum jumps is independent of the positions of the Gea- 
Banacloche states. Exactly as in section IIII A 21 it is 
given by ip( T ) — n^e~ nlT . Within this approximation, 
(|39|l can be computed explicitely even for a time depen- 
dant A9. An elegant way to get the result consists into 
rewriting the sum over the number of quantum jumps 
in a completely different way which does not singularizcs 
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any specific time: 



C. Rabi oscillations in the presence of dissipation 



[] (l + n(0( e j(Ae)(t,) -l))) 



0<t'<l 



(40) 

where n(t') = if no event occurs at time t' and 
n{t') = 1 when a quantum jump occurs at time t' and 
(A9) = 8 m+ — 9 m _. Note that in the above expression, 
t' is not the time of a quantum jump. The formal infi- 
nite product in the r.h.s of (|4H|) can then be expanded 
leading to an expansion involving multi-time correlators 
(n(ii) . ■■n(t' r )) where < t[ < . . . < t' r < t (here, r is 
not the number of quantum jumps). Because n(r) = 
or 1: 



(n(t' 1 )...n(t' r )) = (n(t[))xl[P(t' l+1 \t' l ) 



(41) 



i=i 



where P (t' l+1 \t\) = Prob (n(t' l+1 ) = l\n{t\) = l). At 
short times, quantum jumps provide a renewal process 
and therefore the conditional probabilities P {t'i + i\t'i) are 
directly related to the average density of jumps S(t): 
P{t'i+iK) = S(t[ +1 -t[). Moreover, (n(t[)) = S(t[). 
Here the density of events is a constant: S(t) = jn. The 
final result for the dccohcrcncc functional is thus: 



F[0 m+ ,6m-] = exp 



7fijf ( e j(Ae)(T) - l)dA . (42) 



Finally, expression l|42|l can be interpreted as resulting 
from the accumulation of dccohcrcncc coefficients over in- 
finitesimal periods of time. In the mcsoscopic regime, be- 
cause of the strong coupling regime, the evolution of the 
atoms + cavity system is a forced evolution of the cavity 
state driven by the atomic polarizations \D m (t)), leading 
to the motion of the Gea-Banacloche states \ijj m (t)). For 
each Gea-Banacloche state, this forced motion between 
r and r + dr leaves an imprint in the environment of the 
cavity. The overlap between imprints left by two distinct 
Gea-Banacloche states is precisely the decoherence coeffi- 
cient. During infinitesimal time dr, the imprint left in the 
environment by each Gea-Banacloche state under consid- 
eration is the same as the one left by coherent states with 
the same average number of quanta but a time-dependent 
phase separation (A#)(r). Therefore, the corresponding 
decoherence coefficient is given by: 

D(t, t + dr) ~ cxp (jn (e i(A8)(r) - 1) dr) . (43) 

Since the environment is Markovian, the infinitesimal de- 
coherence coefficients Q43JI associated with different time 
windows [i, t + t] (r» r c ) accumulate through time evo- 
lution, leading to (|42|) . 



1. Spontaneous revivals 

It is now straightforward to compute the Rabi oscilla- 
tion signals by introducing the decoherence coefficients 
for all pairs of Gea-Banacloche states that can appear 
in the reduced atomic density operator. Because we are 
dealing with experimental situations such that the dura- 
tion of experiment is small compared to the cavity dis- 
sipation time, we shall assume that the average number 
of quanta remains equal to n in this paragraph and the 
following ones. 

The result for the Rabi oscillation signal is (q = m + — 
TO_): 



r n 



. (t) Tin 



it) 



(44) 

where V m+i ,„_{t) and lZ m+ ,m-(t) are respectively given 
by equation l|18fl and l|2U|l . Decoherence is contained in 



»-(*) 



-<i q (t)+iQ q (t) w j 1 j c i 1 can 5 e evaluated using 



(|4*2*j) . thus leading to (</> = gt/2\/n): 



d q (t) 
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3/2 



sin , 

gq V 2 



(45) 
(46) 



Note that in these results, n 3 ^ 2 j/g is the dimensionless 
parameter that characterizes the strength of decoherence. 

An ansatz for the upper and lower envelopes of the 
Rabi oscillation signal in the presence of dissipation can 
then be obtained along the lines of section III D 21 

P+(t) = A + ]T|^(i)A|e- d ' (t) (47) 



P_(t) - A + ^2(-iy\1Z q (t)A q \e- d ^ . (48) 



2. Induced revivals 

Rabi oscillation signals in an echo experiment can also 
be computed within the mesoscopic approximation. The 
percussional echo pulse is applied to the atoms at time 
t n . It corresponds to the operator C/ T = ^jL^. Using 
this operator, equation Q can be derived for the case 
of TV atoms with Htc m place of the Jaynes-Cummings 
Hamiltonian. Thus, exactly as for TV = 1, the evolu- 
tion of the atoms + cavity system is reversed after time 
t % . Within the mesoscopic approximation, this means 
that atomic polarization as well as Gea-Banacloche states 
move backward towards their initial positions. The time 
dependence of the associated phases m (T) associated 
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with Gca-Banacloche states is given by: 
J 6 m ( T ) = mgr/2y/E for < r < t v 



1 # m (r) = mg(2t v - r)/2y/n for <r 

Eq. ((42(1 leads to the following decoherence co 

cient (t > ijr) which we write as Fm^m- (t n , t) 
e -d q (t„ ,t)+iQ q {t 7r ,t) w i iere = m+ _ m _) : 

/ 9 sin (arh-\ — sin 



(49) 



(50) 



5 

where ^ = gt 7r /2^/E and 

G q (U,t) = -I (2 sin 2 (?^/2) - sin 2 (g(2^ - 0)/2)) 

^ (51) 
Because of the perfect time reversal for the atoms + cav- 
ity system, the overlap factor lZ m+ ,m_ in the echo exper- 
iment can be expressed in terms of the overlap factor un- 
der free evolution for t > U: n (ccho) (t ni t) = K{2t n - t). 

3. Extension to finite temperature 

It is known that increasing the temperature lowers the 
decoherence time. For the harmonic oscillator, initially in 
a coherent state, the exact solution to the quantum mas- 
ter equation 1(26(1 is well known [23]. At time t, the state is 
no longer pure but appears to be a thermal state with av- 
erage number of quanta n(t) = ny(l — e~ 7 ') translated in 
phase space by ae -7 */ 2 (ht denoting the average number 
of quanta at equilibrium at temperature T). Nevertheless 
for times much shorter than dissipation, thermalization 
can be neglected: the state of the oscillator can still be 
considered as coherent. The analysis of the solution of 
the master equation at finite temperatures shows that 
decoherence at short times is still exponential. The ef- 
fect of temperature is to enhance the damping rate by 
a factor 2n t h + 1 = coth (f3fkjj /2). This suggests that 
the imprint of the superposition of two coherent states 
in the environment during an infinitesimal time inter- 
val dr at finite temperature is obtained by substituting 
7 — > 7 coth ((3Tiu}q/2) in eq. 1)43(1 . Following the previous 
line of reasoning (end of section llll B 2() . the decoherence 
coefficient at time t for a superposition of two coherent 
states is again obtained by summing the decoherence co- 
efficients associated with infinitesimal time intervals be- 
tween r = and t = t. 

This result can be used to derive the evolution of the 
atoms + cavity density matrix at short times and mod- 
erate finite temperature. As long as we can neglect the 
thermalization, the only effect of dissipation is to damp 
the coherences between states \D m± (t)} (g> \ip m± {t)) for 
m + 7^ m„ . As in the zero temperature case, each of the 
state \D m± (t) <g> |^ rra± (t)) is expected to evolve accord- 
ing to the atoms + cavity interaction and the echo pulse 



applied to the system (if any). As explained in the pre- 
vious paragraph, the decoherence coefficient J r m+tTn _(t) 
to be used in eq. 1(44(1 is obtained by replacing 7 by 
7 coth (flhuJo/2) in equations l(45|l . ((46(1 for the free evo- 
lution and ((50(1 . 1(51(1 for the echo experiment. Note that 
this ansatz is expected to be valid only for low tempera- 
tures and at short times such that jt coth (/3huio/2) <C 1. 



IV. DISCUSSION OF THE RESULTS 

A. Method and parameters 

We have considered the Rabi oscillation signal in the 
presence of dissipation for N ~ 1 to TV = 3 atoms, values 
that can be realistically reached in state-of-the-art cav- 
ity QED experiments. Photon numbers n = 10 and 15 
have been considered. All our computations have been 
performed for values of g/j corresponding to the present 
ENS experiment 0. The best published cavity damp- 
ing time is 1 ms (quality factor Q = 3.2 x 10 8 ). Pre- 
liminary tests of an improved experimental setup have 
shown damping rates of 14 ms (Q = 4.5 x 10 9 ) and even 
115 ms (Q = 3.7 x 10 10 ) and these results are to be sub- 
mitted in the near future. Thus, the values of 14 ms 
(3/7 ~ 4310), 5 ms (3/7 ~ 1540 and Q = 1.6 x 10 9 ) 
and 1 ms (g/j ~ 308) for the damping rate of the cavity 
have been considered in our simulations. We focus on the 
case of a zero-temperature bath which can be realistically 
reached as shown in j2^ . The effect of finite temperature 
will be briefly discussed in section llVDI 

Results of the analytical approach described in sec- 
tion llll Ol have been compared to a quantum Monte-Carlo 
simulation of the atoms + cavity system evolution in the 
spirit of [13 . For these simulations, the Adams-Bashford 
scheme of order four has been used to compute the evo- 
lution of the wave function between quantum jumps. 

We first present our results relative to the free evolu- 
tion of the atoms + cavity (spontaneous Rabi oscillations 
revivals) in section TlV Bl and for the echo experiments in 
section IIV CI Consequences of these results for cQED 
and circuit-QED experiments are then discussed in sec- 
tion |IVD2] 



B. Free evolution 

It is interesting to assess the possibility of observing 
spontaneous Rabi oscillation revivals, since, for N > 1, 
such a revival might be observable at shorter times than 
in the N = 1 case. Figure [7| presents a comparison be- 
tween the Rabi oscillation signals resulting from the in- 
teraction with a mesoscopic coherent state containing 15 
photons in average in the dissipationless case and for dis- 
sipation times equals to 1 ms, 5 ms and 14 ms (T = K). 
Figurcs[H]and[!5]present the same comparison for the cases 
of N = 2 and N = 3 atoms respectively. Note that our 
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FIG. 7: Influence of dissipation on the spontaneous revivals 
of the Rabi oscillation signal P = P m =.j (as a function of 
cf> = gt/2y/n) for one atom and m — 1/2 (atom excited) with 
n = 15 photons initially. The graph on the left depicts the 
analytical envelopes P±(t) for (a) no dissipation, (b) 7" 1 ~ 
14 ms, (c) 7 _1 ~ 5 ms and (d) 7 _1 ~ 1 ms. The right part 
of the figure presents the associated Rabi oscillation signals 
obtained from quantum Monte-Carlo simulations (plain lines) 
as well as the associated analytical enveloppes. The vertical 
dashed line corresponds to the largest reachable <j> (atoms at 
100 ms" 1 ). 



analytical model (eqs. (|47|) and (|48|l ) predicts the upper 
and lower envelopes of the Rabi oscillation signal with 
a rather good precision in the presence of dissipation. 
Note that in the N = 1 case, this simulation shows that 
it is not possible to observe spontaneous Rabi oscillation 
revivals with the present cavity (7 < 1 ms). The same 
conclusion is valid for N = 2 and TV = 3: even the partial 
revivals that occur before the main one <f> ~ 2ir should 
not be observable. For 7" 1 = 1 ms, decoherence trans- 
forms the entangled atom + cavity state into a statistical 
mixture before any pair of Gea-Banacloche components 
of the field overlap. As figures [7| to [5] show, improving 
the quality factor of the cavity could enable a direct ob- 
servation of the spontaneous revivals in the presence of 
n = 15 photons. 

However, the experimental apparatus sets a tight lim- 
itation on the interaction time of the atoms with the 
cavity field. The number of atoms flying through the ap- 
paratus at velocities lower than 100 m s _1 is too small in 
a thermal atomic beam to be used in practice. Atomic 
spontaneous emission is another limitation for very slow 
atoms. We stick here to the available apparatus and set 
an upper limit on <f> which could roughly be estimated 
as m ~ 2tt (2.5 /y/n). This upper limit for (f>/2n ranges 
from 0.65 (n = 15) to 0.76 (n = 10) which excludes 
the observation of the complete revival obtained when 
all Gea-Banacloche states overlap again. Nevertheless, 
as can be seen from graphs (b) and (c) on figures [H] and 
13 observing a spontaneous partial revival may be pos- 



FIG. 8: Influence of dissipation on the spontaneous revivals 
of the Rabi oscillation signal P = P m -.j (as a function of 
4> = gt/2\/n) for N = 2 atoms. The initial condition is m = 1 
(all atoms excited) with n — 15 photons. The graph on the 
left depicts the analytical envelopes P±(t) for (a) no dissipa- 
tion, (b) 7 _1 ~ 14 ms, (c) 7 _1 ~ 5 ms and (d) 7 -1 ~ 1 ms. 
The right part of the figure presents the associated Rabi os- 
cillation signals obtained from quantum Monte-Carlo simu- 
lations (plain lines) as well as the associated analytical en- 
veloppes. 



sible for N = 2,3. In the case of two atoms, the signal 
would correspond to the overlap of states \tpi) and IV'-i)- 
In the case of three atoms, the signal would be dominated 
by two overlaps corresponding to IV'1/2) and ^-3/2) 011 

one side and | -0 x / 2 ) and |-0 3 / 2 ) on the other side. Note 

that the partial revival associated with the overlap be- 
tween |f/ ; 3/2) and |^ 3/2) is n °t within reach. Observing 

a spontaneous revival requires <j> ~ tt to remain within 
reach for the slowest atoms. This puts an upper limit on 
the average photon number close to 25 photons. 

Finally, figureEH s h° ws the decay of the modulus of the 
three decoherence coefficients T q (t) as a function of time 
for n = 10 and n — 20 for dissipation times equal to 1, 5 
and 14 ms. It clearly shows that, even if the time of flight 
problem could be circumvented, increasing the number of 
photons puts strong constraints on the dissipation time 
of the cavity. This suggests that working with 10 to 15 
photons in average is a good compromise for observing 
spontaneous partial revivals. 



C. Echo experiment 

In this paragraph, results for the simulation of an echo 
experiment corresponding to an echo pulse at t„ — 30 fxs 
are presented. Figure 1111 shows the simulated echo sig- 
nals for the cQED experiment at LKB obtained with 15 
photons initially and an echo pulse at 30 /xs for N = 1, 
N = 2 and N = 3 atoms and two different values of dis- 
sipation: 14 ms, 5 ms and 1 ms. Note that in the case of 
three atoms, a revival occurs at 150 fis. It corresponds to 
a delayed revival of type (c) on figure^] But although its 




FIG. 9: Influence of dissipation on the spontaneous revivals 
of the Rabi oscillation signal P = P m =.j (as a function of 
4> = gtj2\fn) for N = 3 atoms. The initial condition is m = 
3/2 (all atoms excited) with n — 15 photons. The graph 
on the left depicts the analytical envelopes P±(t) for (a) no 
dissipation, (b) 7 -1 ~ 14 ms, (c) 7 -1 ~ 5 ms and (d) 7 -1 ~ 
1 ms. The right part of the figure presents the associated 
Rabi oscillation signals obtained from quantum Monte-Carlo 
simulations (plain lines) as well as the associated analytical 
enveloppes. 




FIG. 10: Modulus of the decoherence coefficient T q {t) for 
5=1 (plain lines), q — 2 (dashed lines) and q — 3 (dotted 
lines) as a function of <j> for dissipation times 14 ms, 5 ms and 
1 ms in the case of n = 10 photons (left graph) and n = 20 
photons (right graph). The vertical dashed line corresponds 
to the largest reachable cf> (atoms at 100 ms -1 ). 



amplitude makes it visible with a 14 ms dissipation time, 
the time of flight limitation will prevent it from being 
observed in the experiment. The same conclusion holds 
for the delayed revival predicted also in the N = 2 atoms 
case. 



N=l N=2 N=3 
1 1 1 1 1 r 




FIG. 11: Echo signal (function of <f> = gtjlsfn) simulated for 
n = 15 at t-n = 30 fjs and g/2ir — 49 kHz. All graphs on a row 
correspond to the same dissipation time: (a) 7 -1 = 14 ms, (b) 
7 -1 = 5 ms and (c) N = 3 and 7 -1 = 1 ms. All graphs in a 
column correspond to the same numbre of atoms. The arrow 
shows the time of the echo pulse. The vertical dashed line 
corresponds to the largest reachable <f> value (slowest atoms 
at 100 ms~ ). The initial condition is m = 3/2 (all atoms 
excited) . 

D. Finite temperature 

1. Thermalization procedure 

Experiments with Rydberg atoms in high-quality mi- 
crowave cavities are performed at low temperatures (T ~ 
0.8 — 1.4 K). In order to wash out photons resulting 
from thermal leaks, an erasing procedure using auxiliary 
atomic samples is applied Q . Once the erasing procedure 
has been completed, a coherent field is injected inside 
the cavity. Because of the imperfections of the procedure 
and because of the necessary delay t\ ~ 200 [is between 
the erasing sample and the coherent field injection, this 
creates a translated thermal state, partially thermalized. 
This state would then evolve during time T2 ~ 50 [is 
before the experimental atomic sample enters the cavity. 

In order to model this preparation, we have performed 
a quantum Monte-Carlo simulation involving a finite 
temperature environment and an initial thermalization 
period of duration t p . At the beginning of this prepara- 
tion period, a coherent state is injected in the dissipa- 
tive cavity and evolves decoupled from the atoms during 
time t p . Then the coupling to the atoms is turned on 
to model the experiment. A rather pessimistic estimate 
of the average number of thermal photons per mode of 
the reservoir has been used for this simulation (n t h ~ 0.4 
corresponding to T ~ 2 K at 51 GHz). It has been es- 
timated that thermal fluctuations left by the imperfect 
erasure procedure correspond to at most n ~ 0.15 pho- 
tons per mode on the average. Using no = n t h (1 — e -7 * p ), 
this sets "ftp ~ 0.47 used in our numerical thermalization 
protocol. The injected coherent state before thermaliza- 
tion has an amplitude ao = \fn~e ,tv l' 1 in order to take 
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FIG. 12: Rabi oscillation signals obtained from quantui 
Monte-Carlo simulations for one atom using the thermalizE 
tion protocol producing a thermal state with 0.15 photons 
displaced in Fresnel plane by an amplitude corresponding to 
n = 15 photons, for (a) 7 _1 = 14 ms, (b) 5 ms and (c) 
1 ms. Graph (d): comparison between zero temperature sig- 
nals (dashed line) and finite temperature signals (full line) at 
short times for j^ 1 = 5 ms. Graph (e): comparison between 
zero temperature signals (full line) and finite temperature sig- 
nals (dotted line) for 7 _1 = 5 ms. Dashed lines represent the 
upper and lower envelopes P± obtained by taking tempera- 
ture into account by rescaling 7 1— » (1 + 2n t h) 7. 



into account the exponential decay during the prepara- 
tion phase [n denotes the average photon number when 
the experimental atomic sample in injected inside the 
cavity) . 



2. Finite temperature results 

Figure ^] presents the results of these simulations for 
N = 1 atom. The main curve presents the Rabi oscilla- 
tion signals obtained from a quantum Monte-Carlo sim- 
ulation implementing the thermalization procedure de- 
scribed above for 7 -1 = 1, 5 and 14 ms. In order to 
compare it with the analytical model we have used the 
fact that, in the present case, at times short compared to 
the dissipation time, the main effect of finite temperature 
is to speed up decoherence. This is taken into account by 
replacing 7 by 7(1 + 2nth) in the decoherence functional 
(eqs. I|45f) and (|46|) ). This ansatz is used to compute the 
upper and lower envelopes that appear on fig. 1121 

Note that some features are not reproduced by our an- 
alytical ansatz since the overlap factor we use does not 
take into account thermalization of the quasi-coherent 
state \ip m (t)). As shown in graph (d), the initial collapse 
of Rabi oscillations occurs earlier than at zero tempera- 
ture. On the other hand, the envelope of the spontaneous 
revival is well described by our model (see graph (e)). 
This shows that our analytical approach is quite efficient 
in predicting the contrast of spontaneous revivals of the 




FIG. 13: Rabi oscillation signals obtained from quantum 
Monte-Carlo simulations for N = 3 atoms using the ther- 
malization protocol producing a thermal state with 0.15 ther- 
mal photon number displaced in Fresnel plane with n = 15 
for (a) 7~ x = 14 ms, (b) 5 ms and (c) 1 ms. Graph (d): 
comparison between zero temperature signals (dashed line) 
and finite temperature signals (full line) at short times for 
7 _1 = 5 ms. Graph (e): comparison between zero tempera- 
ture signals (full line) and finite temperature signals (dotted 
line) for 7 _1 = 5 ms. Dashed lines represent the upper and 
lower envelopes P± obtained by taking temperature into ac- 
count by rescaling 7 1— > (1 + 2n t h) 7- 

Rabi oscillation signals even at finite temperature. 

Results for the case of three atoms are presented on 
figure^] The same effects as for N — 1 can be observed 
here. Note that thermal fluctuations do reduce the con- 
trast of spontaneous revivals albeit not enough to make 
them unobservable for 7 -1 = 5 and 14 ms. As shown in 
graph (e) , thermal effects reduce the contrast of the main 
secondary revival from 15 to 11 % for 5 ms dissipation 
time. As mentioned above, improving the dissipation 
time reduces the impact of thermalization. 



E. Discussion of the results 

1. Perspectives for the ENS experiment 

Within the context of cQED experiments performed 
at ENS, our results suggests that spontaneous revivals 
of Rabi oscillations could be observed in an improved 
experimental setup with two or three atoms and an initial 
coherent state containing from 10 to 15 photons on the 
average. A convincing test of the generation of three- 
and four-component Schodinger cat states involving 10 
photons and two or three atoms would require to probe 
the phase distribution using an homodyne method 
The splitting of the initial state into N + 1 separated 
phase peaks followed by the recombination of some of 
them at the time of partial spontaneous revivals would 
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provide an experimental proof of the scenario presented 
in the present paper. 

A more detailed insight into the field dynamics could 
be gained by reconstructing the cavity field Wigncr func- 
tion, using the method proposed by and Davidovich and 
already implemented on a one-photon field in p3 | . 



2. Experimental consequences for circuit QED experiments 

Our analysis can be used to discuss the case of circuit- 
QED experiments [l^. Of course, our model does not 
take into account relaxation nor dephasing of the atoms 
themselves, since they are not relevant for the Rydberg 
atom experiments. In the case of circuit QED, the relax- 
ation and decoherence of Josephson qubits must be taken 
into account to obtain a precise model of the dissipativc 
dynamics. In the litterature, the g/-f ratio is around 20 
[25l |26[. Therefore, our model suggests that observing 
the effects discussed in the present paper in circuit QED 
experiments requires an improve of roughly two orders of 
magnitude on this ratio. 

Increasing g/"f would then require either an an in- 
crease of the resonator quality factor or an increase of 
the qubit/cavity coupling. Using classical clctrodynam- 
ics, an upper bound on the value of g can be estimated to 
be of the order of g/u>o — A y / 2a/e r where a is the fine 
structure constant, e r the relative permitivity of the sub- 
strate and A collects geometrical form factors and capac- 
itance ratios determining the microstrip impedance and 
its coupling to the qubit. For resonators in the 1-10 GHz 
frequency range, values of g/2w of the order of 100 MHz 
are quite realistic and would already provide a factor 5 
on .9/7. In present circuit-QED experiments, detection 
is performed by probing the cavity through transmission 
measurements. This requires to lower the quality factor 
from 10 6 [27J to 10 4 which is still in the strong coupling 
regime but not as deeply as for Rydberg atom experi- 
ments. Less dissipative cavities could then be used pro- 
vided one could improve the qubits properties or use an 
alternative measurement method. Recently, new detec- 
tion schemes based on dynamical bifurcation of Joseph- 
son junctions have been developed and provide high con- 
trast, rapid measurement, low back action and absence 
of on-chip dissipation [23, H^, |2jJ. Thus, although an 
important progress is required for circuit-QED devices, 
the rapid improvement of Josephson device technology 
is very encouraging and provides a strong motivation for 
further theoretical studies. 



V. CONCLUSION 

We have studied the resonant interaction of an en- 
semble of N atoms symmctricaly coupled to a resonant 
mcsoscopic field in a cavity. The interaction between the 
atomic ensemble and the cavity produces an entangled 



state with N + 1 components leading to a rich pattern 
for Rabi oscillation revivals that generalizes the ones ob- 
tained in the N = 1 case. In particular, "fractional" 
spontaneous revivals reflecting partial disentanglement of 
the atom + cavity state are expected to occur earlier than 
the first spontaneous revival in the N — 1 case. 

Dissipation in the cavity will lead to decoherence of 
this mcsoscopic entangled state and wc have proposed 
a simple analytical model that enables us to compute 
the Rabi oscillation signals in the presence of dissipation. 
This model provides simple expressions for the sponta- 
neous Rabi oscillations revivals as well as for the ones 
induced in an echo experiment. Analytical results are 
in good agreement with quantum Monte-Carlo simula- 
tions and provide an intuitive view of the evolution of 
the dissipative atoms+cavity system. We have obtained 
an analytical expression for the accumulated decoherence 
of the N + 1-component Schrodinger cat state resulting 
from the atom + field interaction which could be used 
for tests of decoherence. 

We have shown that in a forthcoming generation of 
cQED experiments, spontaneous revivals of Rabi oscilla- 
tions associated with the recombination of a fraction of 
the N + 1 components of the entangled atoms + cavity 
state should be observable. An improvement by a fac- 
tor 10 in the cavity quality factor as well as the use of 
slow atoms (100 m/s) are required. Our analysis also 
suggests that the situation is not so favorable within the 
context of circuit-QED experiments, due to the measure- 
ment limited quality factor. However, rapid progresses 
in this field are extremely encouraging and motivate fur- 
ther theoretical studies taking into account relaxation 
and dephasing of the qubits. The stochastic wave func- 
tion approach could be used to take these dissipative phe- 
nomenon into account at least for the part due to high 
frequency noise. Dissipative dynamics in the presence 
of a strong low frequency 1// noise cannot be accounted 
for within the framework of the Bloch-Redfield equations 
|3lj . Nevertheless, providing a simple analytical model 
for the dissipative dynamics of a combined qubits + cav- 
ity system taking into account all possible sources of dis- 
sipation would is an interesting challenge. 
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APPENDIX A: MESOSCOPIC APPROXIMATION 

In the classical limit, we expect eigenstates of J x to 
remain uncntanglcd with the cavity state. Therefore, let 
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us start from an initial state: 

\<)= E R^, m \J,m")®\a), (Al) 

ra" — — J 

where a = \ffi and R m .m' = (J, m'\e™ JV \J, to). Since we 
are in the mesoscopic regime, most of the weight of the 
state is concentrated in stable subspaces Tt n with n close 
to h. Apart from the term ml = J, the above sum also 
spreads on the lower dimensional stable subspaces. Since 
we expect this contribution to be exponentially weak in 
the mesoscopic domain, we focus the projection of 
l^m) m ©n^™- Because the effective Hamiltonian <(T2|) 
is written in terms of the J x generator, it is useful to 
decompose \^^) on the basis vectors |A™^). We start 
from 



E" 

k,m' 



fe/2 



n-l i 7 (fc-J+m')\ 
"■m,m'\ m' I 



(A2) 



where the sum over k ranges from J — to' > to +00. 
Shifting the index k into p = k + m! — J > and approxi- 
mating y / n k /k\ by \Jh p /pi enables to do the summation 



over ml using \Zm ) 
leads to: 



Em" 



\x 



This finally 



where \D m (t)) and \ip m (t)) are defined respectively by 
cq. <(T5|) and (fHT| . Equation lfTS|) is recovered by expand- 
ing I J, too) over the atomic polarizations at initial time 

|An(0)) (-J < TO < J). 



APPENDIX B: ROTATION MATRICES 

Matrix elements of SU (2) general elements are given 
for example in |32| . The matrix elements R m ,m' = 
( J, to'| exp {iitJ v /2)\ J, to) needed in the present paper 
are: 



i (J-to)!(J-toQ! 
(J + m)!(J + m')\ 



J-max(m,m ) 



E oJ-k 



{-\y- k {2j -k)\ 



k=0 



2 J - k k\{J - to - k)l{J - to' - k)\ 



Starting from excited atoms and looking for the probabil- 
ity of detecting finally all the atoms in the excited state 
involves: 



1 



Rm.J — (—1)' m Rj,m — 7T7 



(2J)! 



2 J V (J — to)!( J + to)! 



E 



IX^) 



(A3) 



When evolved during time t under the effective Hamilto- 
nian (|12(l . this state becomes: 

=e-*E4r e- 89tm ^ |XW) . (A4) 

We now use the matrix to go back to the usual ba- 



sis: IX 



(p)\ 



(p)\ 



In order to rewrite the 



resulting state as a tensor product of an atomic polar- 
ization and a field state, it is necessary to introduce the 
same approximations as before. We first approximate 
y/(p + 1)1 /pi nP ~ 1. The coefficient of | J, J - I) ® \p + 1) 
is then equal to: 



Tf—l -igmt^p+c 



We then replace e % 9 mt V¥+c by e~' l 9mtVp~+l multiplied by 
the phase factor ^^(VFR-Vp+cJ _ This last phase factor 
is then expanded to first order near h which leads to a 
slowly varying phase e i 9mt(i-c)/2Vri _ rp ne resu ^ m g coe f_ 
ficient now only depends on I and p + I. We then define 
k = p + I and extend its summation range from k = 
to +00, introducing an exponentially small error in the 
mesoscopic regime n, > 2J. The resulting vector is now 
our mesoscopic approximation to ^^(t)). It is obtained 
as a tensor product: 

|**(t)) ~ e -*»WS [£)ro(t)) g ^ m{t)) (A5) 



APPENDIX C: DECOHERENCE AND PHASE 
DIFFUSION 

In the jt > 1 case, the decay of the average photon 
number n(t) = ne" 1 ' prevents the sequence of quantum 
jumps from being a renewal process with a stationary 
waiting times probability distribution. Nevertheless, the 
probability that no quantum jump happens between t 
and t + r knowing that one occured at time t is given by: 

no(*,r) = J] (l-jn(t')dt') 

t<t'<t + T 

= exp (~7j£ ^ + r')dr'V (CI) 

The probability distribution for having quantum jumps 
at times £ and f + r is therefore given by 



dU 
dr 



1 n{t + T)e-~ i $° fl{ - t+T ' )dT ' . (C2) 



The probability for having exactly p quantum 
jumps between and t at times < t± < 
< t p < t is equal to 7- , fo,t](ti, • • • ,t p ) = 

ip(o, ti) ip(t 1} t 2 ) ■ . . t/)(tp-i,tp) n (t p , t - t p ): 

V [0A (h,...,t p ) = e-^o «(*') ««' f[ ( 7 n(^)) . (C3) 

fc=i 

Using the exponential relaxation of the mean photon 
number (h(t) = he rt ), this leads to a Poisson distri- 
bution law for the number N[Q, t] of photons emitted be- 
tween and t with average value n(l — e -7 *). Therefore, 



19 



the decoherence coefficient is equal to: 

(e im ^) = exp (n (e w - 1)(1 - e" 7 ")) (C4) 
which is exactly the solution of the master eq. i|26|) . 

APPENDIX D: EFFECT OF DISSIPATION ON 
THE ATOMS + CAVITY 

We consider the dissipative dynamics of the strongly 
coupled atoms + cavity systems in the mesoscopic 
regime. For weak dissipation 7 -C g, we show that a fam- 
ily of generalized Gea-Banacloche states is stable under 
an effective stochastic dynamics naturally arising from 
these approximations. We then consider the evolution 
of the states |^^) (-J < m < J) and show that each 
of them decoheres on a time scale much longer than the 
decay time of the \^m + ) (®m- I coherence for m+ 7^ m_. 
These results validate the simplified analysis presented 
in section ITlI Bl 

1. Effective stochastic dynamics in the mesoscopic 
regime 

The non hcrmitian Hamiltonian that describes the dy- 
namics of the atoms + cavity system between quan- 
tum jumps is Htc ~ ihja i a/2 where Htc is the Tavis- 
Cummings Hamiltonian © . Replacing Htc by the effec- 
tive Hamiltonian 112f> leads to an effective non hermitian 
Hamiltonian over each subspace H. p : 

H = hg^+lj x -i^L(p + J-J'). (Dl) 

The constant term ihJ^f/2 can be discarded in the evolu- 
tion between quantum jumps since it is canceled by the 
normalization of the state vector. It is then useful to 
decompose the atoms + cavity vector over (B p Ti. p as: 

|*(t)>=JV(i)£-7Tl*p(*)> (D2) 

where \^ p (t)} belongs to H p , a(t) denotes some time 
dependent function and Af(t) is the normalization fac- 
tor. Taking a(t) = a(0) e -7 */ 2 enables to absorb the 
p-dependent part in the non-Hcrmitian term of i|Dl|) . As 
a consequence, each vector \^ p (t)) evolves in H. p under 
the non-hermitian Hamiltonian: 

H p = hg^Jp~T~cJ x + i(Hy/2)J z . (D3) 

The time dependence of a reflects the acquisition of in- 
formation arising from the absence of quantum jumps 
on the cavity. Because of the strong coupling between 
the atoms and the cavity, the state of the atoms is also 
altered and this is why dissipation has to be taken into 
account in the evolution of each \^> p (t)) through the non- 
hermitian term in (|D3|) . Note that this induces a change 



in the probabilities for the cavity to release a photon into 
its environment and thus modifies the cavity relaxation. 

It is possible to solve analytically the resulting 
Schrodinger equation for each |\& p (i)) but the resulting 
expressions are very complicated. In particular, tak- 
ing into account the back action of cavity dissipation 
through the atoms on the statistics of photon emission 
leads to very cumbersome expressions. Nevertheless, in 
the strong coupling regime of cQED (7 <C g), neglecting 
the effect of cavity dissipation on the atoms seems to be 
reasonable and turns out to make the problem much more 
tractable. For values of £ such that n{t) = ne~ 7 * S> 1, 
we expect H p to be dominated by its hermitian part. 
Discarding the non hermitian part in l|D3(l means that 
each |*P P (£)) has the same unitary evolution than in the 
dissipationless case. Therefore, within this approxima- 
tion, the atoms + cavity state evolves in the absence of 
quantum jumps between and t as: 

= <^ fi(t)/2 E ^7T ^ l9t ^~ CJ * |*„(0)> 

p=0 V P ' 

(D4) 

where n(t) = \a(t)\ 2 = |a| 2 e -7 *. 

Let us now consider a state of the form l|A3j) and com- 
pute its evolution during a time t in the absence of quan- 
tum jump. Equation (|D4|) leads to: 

(D5) 

Note that, within the mesoscopic approximation, this 
state can be approximated by a factorized state of the 
form (| AB|> taking into account dissipation through the 
exponential decay of n{t). For yt <C 1, eq. (|A5|) provides 
a good approximation to l|D5(l . 

We now consider the effect of quantum jumps on state 
(ID5() . Within the mesoscopic approximation, the action 

of the creation operator on states \Xm^) can be simpli- 
fied: 

a\X$) = y/p + J- m I J, m) ®\p + J - m - 1) 

* v^i^r 1 '). (D6) 

Using this expression, we see that a | 1 S , m(i))n.j. has the 
same form than (|D5|1 . the phase e - l 9 m Wp+ c m front of 
\X$) being replaced by e -^™*VFFT+^ Thig ghowg that 
all states of the form: 

\9 m [a, {*(«)}]) = £ ^ e-^><) \X$) , 

(D7) 

where denotes the set of all phases 4>(p, t) (p > 0) at 
time t, form a stable class along the stochastic evolution. 
The phases <p(p, t) follow a stochastic trajectory consist- 
ing of smooth deterministic evolution periods between 
quantum jumps. The deterministic evolution is ruled by 
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4>(p, t) = mg^/p + c. The quantum jumps correspond to 
discontinuous steps 4>(p, t + ) — <fi(p + l,t). 



2. Evolution of the states |**} 



Studying the evolution of a superposition of states 
of the form (|A3|) leads to consider the evolution 
of a coherence between and ). Af- 



puted in appendix IO we immediately obtain (jt < 1) 



D 



-igmt(^p + +c- ^/p-+c) 



_<>mM 



ter time t, the stochastic evolution of IvP 



x 



X 



p + ,p-{t) — e 

x e fi7/ *(e 4, *+-'- -i>«t (D12 ) 

where m (r) = gm,T fl^ffi and 77 = (p+ — p-)/2n. The 
second factor in 1D12|) is responsible for decoherence of 
the state |^^) of the atoms + cavity system. Because 
of the amplitude e~"^^ 2 a ^ t ) in the atoms + cavity 



produces an ensemble of projectors of the form 



[a,{* + (i)}]>(* m _M*-(t)}]|. 



Denoting by within \P 



Pm + ,m_ (*) the operator obtained by averaging these pro- 
jectors over the measure given by the stochastic trajec- 
tories, we obtain: 



states, the values of p± that contribute to the sum lie 
\J n(t) and therefore \r]\ <C 1 within 



V- 



< 



it) 



-n(t) 



£ 

(p-.p+) 



a(t) 



p-+p+ 



x £L 



^JpAp^\ 

(t)\x^){x% : < 



(D8) 



the mesoscopic regime and for jt < 1. That's why deco- 
herence of a state I A3l can be neglected as in [TTI |. 

Before moving on the m+ 7^ ?n_ case, it is interest- 
ing to see how his results (section III. A) are recovered 
within the present approach. Following we ignore 
the discrete character of p and replace the finite differ- 
ence equation (|D9|I by partial differential equation. Solv- 
ing this equation can easily be done using the character- 
istics method. Starting from the initial condition \ty 



where D p+iP _(t) denotes the average over all this leads to p m ,m{t) = \^m(t)}{^m(t)\ where 
stochastic trajectories of the relative phase factor 
e i(<Mp-. *)-■/■-(?>+>*)) where the phases (j)±(p,t) are rel- 



ative to the states $ ro± . Because the quantum jump 
process is memoryless, these averages obey the following 
set of coupled first order differential equations: 



l*m(*)> 



e _^lj^a{^_ e _ ie{p . t) (D13) 

p=0 



D p+tP _(t) = jn(t)(D. 



p++i,p_+i 



(t) - D p+tP _(t)) 



and 



■n(m-,m+) n (f\ 



(D9) 



where 



^p" +) = 9 (m+ ^/pVTc - m_ ^p~T~c) . (D10) 

Let us first focus on the case m + = m_ = m. Prob- 
ing the {p+,p-) dependence of D p+p _ (t) gives us an in- 
sight of the decoherence of the initial pure state l 1 ?^) of 
the atoms + cavity system. Decoherence arises from the 
quantum jumps that lead to the spreading of 4>(p + ,t) — 
4>(p-,t). Let us estimate the decoherence coefficient be- 



tween \X: 



(p+)\ 



and IX 



(p-)\ 



within the mesoscopic approx- 
1/2-3 the 
is 



imation. Expanding \Jp + c + 1 ~ ^Jp + c 
phase factor associated with a single jump at time t 3 
gmtj /2^/jTjT — gmtj /2^/pZ which is approximately equal 
to gmtj(p- — p + )/4n 3 / 2 for p± close to n (we assume 
jtj < 1 for simplicity). Thus, for m + = m_ = m, the 
decoherence factor D p+tP _(t) can be approximated by 



D 



p+,p- 



(t) ~ e ~ i 9 mt (.VP+ + c ~ Vp- + c ) / e ' 



2ft ^ i 



<3 2v« ) 

. (Dll) 

Remembering that the statistics of occurrence times of 
quantum jumps for all states of the form l|D7|) is corn- 



er, i)=gm/ ^c + p + n(r) - n(t) dr . (D14) 
Jo 

Evaluating the integral and for all p leads to expressions 
corresponding to eqs. (18. a) to (18. c) of fill ]. 

The case m + 7^ rn_ can then be studied along the same 
lines. Within the mesoscopic approximation, the phase 
jump associated with a quantum jump occuring at time tj 
can be evaluated as exp (igtj(m„ — m + ) /2^Jn{tj)). Note 
that it does not vanish for p + = p_ ~ h(tj). This is why, 
decoherence for m + 7^ m_ occurs on a much shorter time 
scale than the decoherence of the state I n princi- 

ple, p m+ .m^{t) could be computed from the formalism 
presented here but this is not necessary for the present 
purpose. In the jt 1 case, the problem can be sim- 
plified by considering that the evolution of the atoms + 
cavity state with initial condition \^>^) produces a pure 
state as in the previous paragraph and by approximating 
this pure state by ( scc CC L- J35H . The argumen- 

tation presented in section IlII Bl then leads to the deco- 
herence properties of the atoms + cavity state suitable 
for this regime. 
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